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ABSTRACT 


An  effective  computational  strategy  has  been  developed  for  the  analysis  of  large  and 
complex  structures.  The  strategy  is  based  on  generating  the  response  of  the  complex  structure 
using  large  perturbations  from  that  of  a  lower-order  (simpler)  model  associated  with  a  simpler 
structure  (or  a  simpler  mathematical/discrete  model  of  the  original  structure).  The  three  key 
elements  of  the  strategy  are:  (a)  mixed  (or  primitive  variable)  formulation  with  the  fundamental 
unknowns  consisting  of  generalized  displacements  and  stress  parameters;  (b)  operator  splitting, 
or  a  reduction  method  to  relate  the  arrays  and  degrees  of  freedom  of  the  original  complex  struc¬ 
ture  to  those  of  the  simpler  system;  and  (c)  efficient  iterative  process  for  the  generation  of  the 
response  of  the  complex  structure  starting  from  that  of  the  simpler  system.  The  strategy  has  been 
successfully  applied  to  a  number  of  linear  and  nonlinear  stress  analysis  problems,  free  vibration 
and  nonlinear  structural  dynamics  problems.  The  strategy  was  also  used  to  obtain  accurate 
transverse  stresses  in  laminated  composite  plates  and  shells,  using  the  two-dimensional  first- 
order  shell  (plate)  theory  as  the  simpler  model. 

INTRODUCTION 

The  work  under  this  grant  focused  on  the  development  of  an  effective  computational 
strategy  for  the  solution  of  large-scale  structural  problems.  The  strategy  developed  combines  the 
following  three  major  characteristics: 

1)  gives  physical  insight  about  the  response; 

2)  helps  in  assessing  the  adequacy  of  the  computational  model;  and 

3)  is  highly  efficient. 

The  first  characteristic  is  a  direct  consequence  of  the  basic  idea  of  the  strategy;  namely, 
generating  the  response  of  a  complex  structure  using  large  perturbations  from  that  of  a  simpler 
structure  (or  a  sequence  of  simpler  structures).  Examples  of  complex  and  corresponding  simpler 
structural  systems  are  given  in  Table  1. 

The  second  characteristic  is  significant  because  of  the  importance  of  assessing  the 
reliability  of  response  predictions  (i.e.,  the  agreement  between  the  response  predictions  of  the 
computational  model  and  those  of  the  actual  structure).  In  the  absence  of  reliable  and  practical 
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error  measures,  information  about  the  sensitivity  of  the  response  to  modeling  details  can  help  in 
identifying  the  regions  of  questionable  accuracy,  and  in  the  adaptive  refinement  of  the  comput¬ 
ational  model.  Sensitivity  information  can  be  obtained  by  using  more  than  one  mathematical 
and/or  discrete  model  of  the  structure.  Only  the  lower-order  model  equations  are  solved.  The 
higher-order  model  arrays  are  used  to  generate  the  sensitivity  derivatives. 

High  computational  efficiency  is  achieved  by  reducing  the  time  required  for  the  accurate 
numerical  simulation  of  the  response  which,  in  turn,  requires: 

a)  reducing  the  number  of  degrees  of  freedom  used  in  the  initial  discretization  (the  initial 
discrete  model  is,  in  many  cases,  dictated  by  the  topology  rather  than  the  complexity  of  the 
response); 

b)  exploiting  the  major  features  of  the  new  and  emerging  computing  systems  (viz.,  vector, 
parallel,  and  AI  capabilities).  Al-knowledge  based  systems  are  used  in  the  initial  selection  and 
adaptive  refinement  of  the  model  as  well  as  in  the  selection  of  numerical  algorithms.  Work  in  the 
last  period  focused  only  on  reducing  the  number  of  degrees  of  freedom  used  in  the  initial  dis¬ 
cretization. 

The  three  key  elements  of  the  strategy  are  as  follows:  1)  use  of  mixed  (or  primitive 
variable)  formulation  with  the  fundamental  unknowns  consisting  of  generalized  displacements 
and  stress  parameters;  2)  a  novel  operator  splitting  technique  which  allows  the  restructuring  the 
governing  equations  of  the  complex  structure  in  such  a  manner  as  to  delineate  the  contributions 
to  the  simpler  system;  3)  application  of  an  iterative  procedure  or  a  reduction  method  for  the 
efficient  generation  of  the  response  of  the  complex  structure  starting  from  that  of  the  simpler 
system. 

Two  general  approaches  have  been  developed  for  selecting  the  simpler  model  and  estab¬ 
lishing  the  relations  between  the  original  and  simpler  models,  namely,  hierarchical  modeling, 
and  decomposition  or  partitioning  strategy.  Successful  applications  of  the  strategy  to  static,  free 
vibration  and  dynamic  response  problems  of  composite  panels  have  been  reported  in  the  publica¬ 
tions  listed  subsequently.  Also,  application  has  been  made  of  the  strategy  to  three-dimensional 
analysis  of  multilayered  composite  panels.  Both  stress  and  free  vibration  problems  have  been 
considered. 
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BASIC  IDEA  AND  KEY  ELEMENTS  OF  STRATEGY 
Application  of  the  strategy  developed  to  the  static  stress  and  vibration  analyses  of  both 
symmetric  and  unsymmetric  structures  has  been  reported  in  Refs.  1  to  7.  Herein,  the  application 
to  nonlinear  static  problems  is  outlined.  For  convenience,  simpler  structure  is  chosen  as  a 
structure  with  symmetric  geometry  and  symmetric  (orthotropic)  material  characteristics.  The 
unsymmetric  nonlinear  response  of  the  original  anisotropic  structure  is  approximated  by  a  linear 
combination  of  symmetric  and  antisymmetric  response  vectors  (which  are  uncoupled  in  the 
simpler  structure).  Each  of  the  symmetric  and  antisymmetric  vectors  is  obtained  by  using  only  a 
fraction  of  the  degrees  of  freedom  of  the  original  Finite  element  model,  as  briefly  described 
subsequently. 

1.  Governing  Finite  Element  Equations 

A  total  Lagrangian  formulation  is  used,  and  the  structure  is  discretized  by  using  two-field 
mixed  finite  element  models.  The  loading  is  controlled  by  a  single  parameter  q.  The  governing 
finite  element  equations  describing  the  large  deflection,  nonlinear  response  of  the  structure  can 
be  written  in  the  following  compact  form: 


J/h (H,X,q)l  r-|F|  IShjHl  +  J  G’(X)  1  |.l 
|/x(H,X,q)J  LlS|*  .  JW  +W(H,X)J  l|[pJ 


(1) 


where  (H)  is  the  vector  of  stress-resultant  parameters;  {X}  is  the  vector  of  nodal  displacements; 
|F|  is  the  global  flexibility  matrix;  |S|  is  the  strain-displacement  matrix;  { P }  is  the  normalized 
load  vector;  |G(X)}  and  { Af(H,X) }  are  vectors  of  the  nonlinear  terms  for  the  structure;  super¬ 
script  t  denotes  transposition;  and  a  dot  (.)  refers  to  a  zero  submatrix,  or  a  zero  subveclor. 

The  determination  of  the  deformed  configuration  of  the  structure  corresponding  to  speci¬ 
fied  values  of  the  parameter  q  is  accomplished  by  solving  the  nonlinear  system  of  algebraic 
equations,  equations  (1),  using  an  incremental-iterative  technique  such  as  the  Newton-Raphson 
method.  The  recursion  formula  for  the  rth  iteration,  can  be  written  in  the  following  form: 

LlS|'+[s"»]'  [k'"]  _  W  1/xJ 
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and  the  range  of  i  is  1  to  the  total  number  of  stress-resultant  parameters  in  the  model. 

In  order  to  reduce  the  cost  of  generating  the  deformed  configurations  of  the  structure, 
reduction  methods  have  been  developed  for  substantially  reducing  the  number  of  degrees  of 
freedom  used  in  the  initial  discretization.  These  methods  are  based  on  successive  applications  of 
the  finite  element  method  and  the  classical  Bubnov-Galerkin  technique.  The  finite  element 
method  is  used  to  generaie  a  few  global  approximation  vectors  (or  modes).  The  Bubnov- 
Galerkin  technique  is  then  used  to  approximate  the  nonlinear  equations,  equations  (1),  by  a  small 
set  of  nonlinear  equations  in  the  amplitudes  of  these  modes.  An  effective  set  of  modes  was 
found  to  be  the  path  derivatives  (i.e.,  the  various-order  derivatives  of  the  response  quantities  with 
respect  to  a  control  parameter  (such  as  load,  displacement,  or  arc-length  parameter  in  the  solution 
space).  The  equations  used  in  evaluating  the  path  derivatives  are  obtained  by  successive  dif¬ 
ferentiation  of  the  original  nonlinear  equations,  equations  (1),  with  respect  to  the  control 
parameter.  The  left-hand  side  matrix  in  these  equations  is  the  same  as  that  of  equations  (2).  The 
details  of  application  of  reduction  methods  to  the  generation  of  the  equilibrium  path  are  given  in 
Refs.  8,  9  and  10,  and  involve  the  following  basic  steps: 

i)  Evaluation  of  the  path  derivatives,  and  formation  of  the  reduced  equations  which  ap¬ 
proximate  the  original  nonlinear  equations 

ii)  Generation  of  the  approximate  solutions  associated  with  different  values  of  the  control 
parameter,  using  an  incremental/iterative  approach,  in  conjunction  with  the  reduced  equations. 

iii)  Sensing  the  error  resulting  from  the  use  of  the  approximate  reduced  equations.  When¬ 
ever  the  error  exceeds  the  prescribed  tolerance,  an  iterative  procedure  is  used  to  generate  an 
improved  solution,  a  new  (updated)  set  of  path  derivatives,  and  a  new  set  of  reduced  equations. 

The  most  time-consuming  steps  of  the  procedure  are  those  associated  with  visiting  the 
original,  full  system  of  equations,  namely,  generation  of  an  initial  (or  improved)  nonlinear 
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solution,  and  evaluation  of  path  derivatives.  Hence,  the  reduction  in  the  size  of  the  analysis 
model  used  in  these  steps,  as  outlined  in  the  succeeding  subsection,  can  reduce  the  total  comput¬ 
ational  effort. 

2. 

To  simplify  the  nonlinear  analysis  of  anisotropic  structures,  each  of  the  response  vectors  in 
equations  (1)  and  (2)  is  decomposed  into  symmetric  and  antisymmetric  components  (with  the 
same  type  of  symmetry  as  that  exhibited  by  the  structure;  or  as  desired  in  order  to  reduce  the  size 
of  the  analysis  model). 

Each  of  the  arrays  in  equations  (1)  and  (2)  is  partitioned  into  the  sum  of  the  array  associ¬ 
ated  with  the  corresponding  symmetric  response  plus  correction  terms.  The  resulting  equations 
are  embedded  into  a  single  parameter  family  of  equations  of  the  form  given  subsequently: 


The  recursion  formula  for  the  rth  iteration,  equations  (4),  can  be  cast  in  the  following  form: 

fr  -ifi0  isu[s<«]„i  r  ~iFix  isil+[s«]x-i\fAH,« 

vjjsi:,+[s('l  Ik"1],,  J  Lisix+Is"’]'*,  Ik"’],.  J/axI 


l/xl„  IfxJx 


(6) 


In  equations  (5)  and  (6),  X  is  a  tracing  parameter  which  identifies  all  the  correction  terms 
(temis  with  subscript  X)  to  the  arrays  associated  with  the  symmetric  response  (terms  with  sub¬ 
script  o).  For  X=(),  the  structure  has  a  symmetric  response,  and  for  X=l,  the  response  is 
unsymmetric. 

The  equations  used  in  evaluating  the  path  derivatives  can  also  be  embedded  in  a  single* 
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parameter  family  of  equations  and  cast  in  a  form  similar  to  that  of  equations  (6). 

3.  Generation  of  the  Response  Vectors 

The  solutions  of  equations  (5)  and  (6)  at  A.=0  are  used  in  conjunction  with  either  a)  a 
reduction  method,  or  b)  preconditioned  conjugate  gradient  (PCG)  technique  to  generate  the 
corresponding  solutions  at  A=l.  This  results  in  reducing  the  size  of  the  analysis  model  required 
for  the  generation  of  these  solutions.  A  detailed  description  of  the  PCG  technique  is  given  in 
Refs.  2  and  6.  In  this  section  a  brief  outline  of  the  reduction  method  is  given. 

The  method  is  based  on  generating  a  few  global  approximation  vectors  at  X=().  The 
solutions  of  equations  (5)  and  (6)  are  then  expressed  as  linear  combinations  of  these  global 
approximation  vectors  in  the  following  transformation: 


where  (Th)  and  |  Tx  |  are  transformation  matrices  and  |\j/j  is  a  vector  of  undetermined  coeffi¬ 
cients  (amplitudes  of  global  approximation  vectors)  which  are  functions  of  X. 

An  effective  choice  for  the  global  approximation  vectors  (the  columns  of  the  matrices  |P|) 
consists  of  the  solution  corresponding  to  X=()  and  its  various-order  derivatives  with  respect  to  X 
(evaluated  at  A,=0),  i.e.. 


The  total  number  of  approximation  vectors  is  r,  which  is  considerably  smaller  than  the  total 
number  of  degrees  of  freedom  in  the  model. 

The  equations  used  in  generating  the  global  approximation  vectors  are  obtained  by  succes¬ 
sive  differentiation  of  the  original  equations  (Equations  (3)  and  (6))  with  respect  to  X.  For 
symmetric  loading,  the  vectors  and{X}x=o  and  all  their  even-order  derivatives  with 

respect  to  X  are  symmetric,  and  all  the  odd-order  derivatives  are  antisymmetric.  The  Bubnov- 
Galerkin  technique  is  then  used  to  approximate  the  original  equations  (equations  (5)  or  (6))  by  a 
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reduced  system  of  equations  in  the  unknown  parameters  (components  of  the  vector  {y}).  The 
details  of  application  of  the  strategy  to  eigenvalue  problems  is  described  in  Refs.  3  and  7.  Its 
application  to  nonlinear,  large  deflection  problems  is  described  in  Refs.  I  and  6. 

IMPORTANCE  OF  THE  COMPLETED  STUDY 
Because  of  the  requirements  of  high  performance,  light  weight  and  economy  for  future 
flight  vehicles,  and  the  associated  stringent  design  criteria,  the  prediction  of  the  response  of  these 
structures  is  likely  to  require  more  sophisticated  analysis  models  than  has  heretofore  been  done. 
Also,  analysis  is  needed  to  reduce  the  dependence  on  extensive  and  expensive  testing  which  is 
frequently  component  or  mission  oriented.  Therefore,  there  is  a  need  for  innovative  solution 
methodologies  and  effective  computational  strategies,  such  as  the  one  developed  under  the 
present  grant,  for  handling  large-scale  structural  and  coupled  problems.  The  work  completed 
under  the  present  grant  is  a  first  step  towards  filling  this  need  and  coupling  the  physics  of  the 
problem  with  the  solution  strategy. 
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Table  1  -  Examples  of  Complex  and  Corresponding 
Simpler  Structural  Systems 


Original  Structure 

Simpler  Structure 

Structure  with  complicated 
geometry 

•  Structure  with  simpler  geometry 
(see  Fig.  1) 

Anisotropic  structure 

♦  Orthotropic  structure 

Symmetric  structure  with 
unsymmetric  boundary 
conditions 

•  Symmetrized  structure  with 
SYMMETRIZED  boundary 
conditions 

Three-dimensional  elasticity 
model  of  laminated  shell 

•  Two-dimensional  first-order 
shear  deformation  model 

Stiffened  shell 

•  Unstiffened  shell 

Fine  (enriched)  grid  model 

•  Coarse  grid  model 
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Computational  Strategy  for  Analysis  of  Quasi- 
Symmetric  Structures 

By  Ahmed  K.  Noor,1  Member,  ASCE,  sod  Samira  L.  Whitworth1 

As* tract:  An  efficient  computational  strategy  is  presented  for  reducing 
the  cost  of  generating  the  response  of  quasi-symmetric  structures.  The 
three  key  elements  of  the  strategy  are:  a)  Use  of  mixed  finite  element 
models  having  independent  shape  functions  for  the  internal  forces 
(stress  resultants)  and  generalized  displacements  with  the  internal  forces 
allowed  to  be  discontinuous  at  inierciemcnt  boundaries;  b)  operator 
splitting,  or  additive  decomposition  of  the  different  arrays  in  the 
governing  equations  into  the  contributions  to  a  symmetrized  response 
plus  correction  terms;  and  c)  application  of  a  preconditioned  conjugate 
gradient  technique  to  generate  the  unsymmetric  response  as  the  sum  of 
symmetric  and  antisymmetric  modes,  each  obtained  using  approximately 
half  the  degrees  of  freedom  of  the  original  model.  The  preconditioning 
matrix  is  taken  to  be  the  matrix  associated  with  the  symmetrized 
response.  The  effectiveness  of  the  proposed  strategy  is  demonstrated  by 
means  of  two  numerical  examples  of  an  anisotropic  shallow  panel  with 
a  quadrilateral  planform,  and  an  anisotropic  conical  panel.  Also,  the 
potential  of  the  proposed  strategy  for  solving  nonlinear  problems  of 
quasi-symmetric  structures  is  discussed. 


'Prof  ,  George  Washington  Univ..  Mail  Stop  269,  NASA  I  j,nglry 
Hampton.  VA  23665 
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'Programmer  Analyst,  George  Washington  Univ.,  NASA  Langley  Res  Ctr 
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Note.  Discussion  open  until  August  1,  1988.  To  extend  the  closing  date  one 
month,  a  written  request  must  be  filed  with  the  ASCE  Manager  of  Journals  The 
manuscript  for  this  paper  was  submitted  for  review  and  possible  publication  on 
Pebruary  1 1 ,1987.  This  paper  is  part  of  the  JotmmJ  of  Erngimotrirng  Mtckanici  Vol 
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Model-size  reduction  for 
the  non-linear  dynamic 
analysis  of  quasi- 
symmetric  structures 


Ahmed  K.  Noor  and  Jeanne  M.  Peters 

George  Washington  University.  NASA  Langley 
Research  Center.  Mail  Stop  269,  Hampton, 
Virginia  23665.  USA 
(Received  March  1987) 


ABSTRACT 

A  computational  procedure  is  presented  for  the  efficient 
non-linear  dynamic  analysis  of  quasi-symmetric 
structures.  The  procedure  is  based  on  approximating  the 
imsymmetric  response  vectors ,  at  each  time  step,  by  a  linear 
combination  of  symmetric  and  antisymmetric  vectors,  each 
obtained  using  approximately  half  the  degrees  of  freedom 
of  the  original  model.  A  mixed  formulation  is  used  with 
the  fundamental  unknowns  consisting  of  the  internal 
forces  (stress  resultants),  generalized  displacements  and 
velocity  components.  The  spatial  discretization  is  done  by 
using  the  finite  element  method,  and  the  governing  semi- 
discrete  finite  element  equations  are  cast  in  the  form  of 
first-order  non-linear  ordinary  differential  equations.  The 
temporal  integration  is  performed  by  using  implicit 
multistep  integration  operators.  The  resulting  non-linear 
algebraic  equations,  at  each  time  step,  are  solved  by  using 
iterative  techniques.  The  three  key  elements  of  the 
proposed  procedure  are:  (a)  use  of  mixed  finite  element 
models  with  independent  shape  functions  for  the  stress 
resultants,  generalized  displacements,  and  velocity 
components  and  with  the  stress  resultants  allowed  to  be 
discontinuous  at  interelement  boundaries;  (b)  operator 
splitting,  or  restructuring  of  the  governing  discrete 
equations  of  the  structure  to  delineate  the  contributions 
to  the  symmetric  and  antisymmetric  vectors  constituting 
the  response;  and  (c)  use  of  a  two-level  iterative  process 
(with  nested  iteration  loops)  to  generate  the  symmetric 
and  antisymmetric  components  of  the  response  vectors  at 
each  time  step.  The  top-  and  bottom-level  iterations 
(outer  and  inner  iterative  loops)  are  performed  by  using 
the  Newton  Raphson  and  the  preconditioned  conjugate 
gradient  (PCCi)  techniques,  respectively.  The 
effectiveness  of  the  proposed  strategy  is  demonstrated  by 
means  of  a  numerical  example  and  the  potential  of  the 
strategy  for  solving  more  complex  non-linear  problems  is 
discussed. 
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NONLINEAR  DYNAMIC  ANALYSIS  OF  QUASI-SYMMETRIC 
ANISOTROPIC  STRUCTURES 

Ahmed  K.  NooRf  and  Jeanne  M.  Peters! 

George  Washington  University  at  NASA  Langley  Research  Center,  Hampton,  VA  23665,  U.S.A. 

Abstract — An  efficient  computational  strategy  is  presented  for  the  nonlinear  dynamic  analysis  or 
quasi-symmetric  anisotropic  structures.  A  mixed  formulation  is  used  with  the  fundamental  unknowns 
consisting  of  stress  resultants,  generalized  displacements  and  velocity  components.  The  governing 
semi-discrete  finite  element  equations  consist  of  a  mixed  system  of  algebraic  and  ordinary  differential 
equations.  The  temporal  integration  of  the  differential  equations  is  performed  by  using  an  explicit 
half-station  central  difference  method.  The  three  key  elements  of  the  strategy  are:  (a)  use  of  mixed  finin' 
element  models  with  independent  shape  functions  for  the  stress  resultants,  generalized  displacements  and 
velocity  components  and  with  the  stress  resultants  allowed  to  be  discontinuous  at  interelement  boundaries; 
(b)  operator  splitting,  or  additive  decomposition  of  the  different  arrays  in  the  governing  equations  into 
the  contributions  to  a  symmetrized  response  plus  correction  terms;  and  (c)  application  of  a  preconditioned 
conjugate  gradient  technique  to  generate  the  unsymmetric  response  of  the  structure,  at  each  time  step, 
as  the  sum  of  symmetric  and  antisymmetric  modes,  each  obtained  using  approximately  half  the  degrees 
of  freedom  of  the  original  model.  The  preconditioning  matrix  is  taken  to  be  the  matrix  associated  with 
the  symmetrized  response. 

The  effectiveness  of  the  proposed  strategy  is  demonstrated  by  means  of  a  numerical  example  and  the 
potential  of  the  proposed  strategy  for  solving  more  complex  nonlinear  problems  is  discussed. 
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Parallel  Processing  in  Finite  Element  Structural  Analysis 


Ahmed  K.  Noor 

Cieorge  Washinglon  University,  NASA  Langley  Research  Center,  Hampton,  Virginia 


Abstract.  A  brief  review  is  made  of  the  fundamental  concepts 
and  basic  issues  of  parallel  processing.  Discussion  focuses  on 
mechanisms  for  parallel  processing,  construction  and  implemen¬ 
tation  of  parallel  numerical  algorithms,  performance  evaluation 
of  parallel  processing  machines  and  numerical  algor ilhnts,  and 
parallelism  in  futile  element  compulations.  A  novel  partitioning 
strategy  is  outlined  for  maximizing  (he  degree  of  parallelism  on 
computers  with  a  small  number  of  powerful  processors. 
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A  COMPUTATIONAL  STRATEGY  FOR  MAKING  COMPLICATED 
STRUCTURAL  PROBLEMS  SIMPLE 

Ahmed  K.  NOOR  and  Jeanne  M.  PETERS 

George  Washington  University ,  NASA  Langley  Research  Center,  Hampton,  VA,  U.S.A. 


Received  17  February  1988 
Revised  manuscript  reccivcil  22  June  1988 


An  effective  computational  strategy  is  presented  for  the  analysis  of  large  and  complex  structures. 
The  strategy  is  based  on  generating  the  response  of  the  complex  structure  using  large  perturbations 
from  that  of  a  lower-order  (simpler)  model  associated  with  a  simpler  structure  (or  a  simpler 
mathcmaticnl/discrctc  model  of  the  original  structure).  The  three  key  elements  of  the  strategy  arc:  (a) 
mixed  (or  primitive  variable)  formulation  with  the  fundamental  unknowns  consisting  of  generalized 
displacements  and  stress  parameters;  (h)  operator  splitting,  or  a  reduction  method  to  lelnic  the  arrays 
and  degrees  of  freedom  of  the  original  complex  structure  to  those  of  the  simpler  system;  and  (c) 
efficient  iterative  process  for  the  generation  of  the  response  of  the  complex  structure  starting  from  that 
of  the  simpler  system.  The  effectiveness  of  the  proposed  strategy  is  demonstrated  by  means  of  two 
numerical  examples. 
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Assessment  of  shear  deformation  theories  for 
multilayered  composite  plates 

Ahmed  K  N&or  and  W  Scott  Burton 

George  Washington  University,  NASA  l/uigtey  Research  Center,  Hampton,  VA  23665 

A  review  is  made  of  ihe  different  approaches  used  for  modeling  multilayered 
composite  plates.  Discussion  focuses  on  different  approaches  for  developing 
two-dimensional  shear  deformation  theories;  classification  of  two-dimensional 
theories  based  on  introducing  plausible  displacement,  strain  and/or  stress  as¬ 
sumptions  in  the  thickness  direction;  and  first-order  shear  deformation  theories 
based  on  linear  displacement  assumptions  in  the  thickness  coordinate.  Extensive 
numerical  results  are  presented  showing  the  effects  of  variation  in  the  lamination 
and  geometric  parameters  of  simply  supported  composite  plates  on  Ihe  accuracy 
of  the  static  and  vibrational  responses  predicted  by  six  different  modeling 
approaches  (based  on  two-dimensional  shear  deformation  theories).  The  standard 
of  comparison  is  taken  to  be  the  exact  three-dimensional  elasticity  solutions. 
Some  of  the  future  directions  for  research  on  the  modeling  of  multilayered 
composite  plates  are  outlined. 
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Stress  and  Free  Vibration  Analyses  of  Multilayered 
Composite  Plates 


Ahmed  K.  Noor&  W.  Scott  Burton 

(ieorge  Washington  l  Inivetsily.  NASA  I  angles-  Research  <  enter.  Mail  Stop  2(>'>. 
I  lampion.  Virginia  2.1(i(i5-5225.  I  ISA 


A  IIS  I  II A(  I 

A  two-phase  computational  procedure  is  presented  for  the  accurate  predic¬ 
tion  i>l  the  vibration  frequencies,  stresses  and  deformations  in  multilayered 
composite  plates.  In  the  first  phase  a  two-dimensional  first-order  shear 
tleformalion  theory  is  used  to  predict  the  global  response  characteristics 
(vibration  frequencies,  '<ivcmge’  thmuf’li-the-iliickness  displacements  and 
rotations }  us  well  as  the  in-plane  stress  ami  strain  components  in  the  different 
layers.  In  the  second  phase,  equilibrium  equations  and  constitutive  relations 
of  the  three-dimensional  theory  oj  elasticity  are  used  to:  (I)  calculate  the 
transverse  stresses  and  strains  as  well  as  the  transverse  strain  energy  densities 
in  the  different  layers;  (2)  provide  belter  estimates  for  the  composite  shear 
correction  factors;  and  (  U  calculate  corrected  values  for  the  vibration  fre¬ 
quencies,  displacements,  and  in  plane  strains  anti  stresses 

f  or  simply  supported  plates  the  predictions  of  the  proposal  procedure 
are  shown  to  be  in  close  agreement  with  exact  three dimensional  elasticity 
solutions  for  a  wide  range  of  lamination  and  geometric  parameters.  A  Iso.  the 
potential  of  the  proposed  procedure  for  use  in  conjunction  with  large-st  ale 
finite  element  models  of  composite  strut  lures  is  discussed. 
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ANALYSIS  OH  MULTILAYERED  ANISOTROPIC  PLATES  - 
A  NEW  LOOK  AT  AN  OLD  PROBLEM 


S9-12&4CP 


Ahmed  K.  Noor*  und  W  Scon  Burton  f 
George  Washington  University 
NASA  Langley  Research  Center 
I  lampion,  Virginia 


Abstract 


A  study  is  made  of  the  effects  of  variation  in  the  lamina¬ 
tion  and  geometric  parameters  of  multilayered  anisotropic 
(nononhotropic)  plates  on  the  accuracy  of  the  static  and 
vibrational  responses  predicted  by  six  modeling  approaches, 
bused  on  two-dimensional  shear  deformation  theories.  Two 
key  elements  distinguish  the  present  study  from  previous 
studies  reported  in  the  literature:  a)  the  standard  of  comparison 
is  taken  to  be  the  exact  three-dimensional  elasticity  solutions. 
and  b)  quantities  compared  are  not  limited  to  gross  response 
characteristics  (eg.,  vibration  frequencies,  strain  eitcigy 
components,  average  through-ihe-thickncss  displacements  und 
rotations),  but  include  detailed,  through-lhe-lhickncss,  distribu¬ 
tions  of  displacements,  stresses  and  strain  energy  densities. 

The  modeling  approaches  considered  include  first  order 
shear  deformation  theory  (with  five  displacement  parameters  to 
characterize  the  deformation  in  the  thickness  direction); 
first  order  theory  with  the  transverse  normal  stresses  and 
strains  included  (six  displacement  parameters);  two  higher 
order  theories  (with  II  and  IK  displacement  parameters);  a 
simplified  higher-order  theory  (with  S  displacement 
parameters);  and  a  predictor-corrector  approach,  used  in 
conjunction  with  (he  first-order  shear  deformation  theory  (with 
five  displacement  parameter  in  the  predictor  phase) 

Based  on  the  numerical  studies  conducted,  the  predictor - 
corrector  upproach  appears  to  be  the  most  effective  among  the 
six  modeling  approaches  considered.  For  antisymmetrically 
luminuted  rectangular  plates  the  response  quantities  obtained 
by  the  predictor-corrector  approach  are  shown  to  l>e  in  close 
agreement  with  exact  three-dimensional  elasticity  solutions  lor 
a  wide  range  of  lamination  and  geometric  parameters.  The 
potential  of  this  approach  for  predicting  the  response  of  multi¬ 
layered  anisotropic  plates  with  complicated  geometry  is  also 
discussed. 
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Assessment  of  Computational  Models  for  Multilayered 
Anisotropic  Plates 


Ahmed  K.  Noor  &  W.  Scot!  Burton 

George  Washington  University.  MS-269.  NASA  l.anglcy  Research  ('enter,  Hampton, 

Virginia  2.W»5  5225.  USA 


ABSTRACT 

A  study  is  made  of  the  effects  of  variation  in  llte  lamination  and  geometric 
parameters  of  multilayered  anisotropic  ( nonorthotropic )  plates  on  the 
accuracy  of  the  static  and  vibrational  t  espouses  predicted  by  eight 
modeling  approaches,  hosed  on  two-dimensional  shear  deformation 
theories.  Two  key  elements  distinguish  the  present  study  from  previous 
studies  reported  in  llte  literature  (!)  the  standard  of  comparison  is  taken  to 
be  the  exact  three-dimensional  elasticity  solutions,  and  (2)  quantities 
compared  arc  not  limited  to  gross  response  characteristics  (e  g.  vibration 
frequencies,  strain  energy  components,  average  through  the  thickness 
displacements  and  rotations ),  but  include  detailed  Ihrongli-the-thickness 
distributions  of  displacements,  stresses  and  strain  energy  densities. 

The  modeling  approaches  considered  include  first  order  shear- 
deformation  theory  (with  five  displacement  parameters  to  characterize  the 
deformation  in  the  thickness  direction);  first-order  theory  with  the 
transverse  normal  stresses  and  strains  included  (six  displacement  para¬ 
meters);  two  higher  order  theories  (with  II  and  18  displacement  para¬ 
meters);  a  simplified  higher  order  theory  ( with  five  displacement  para¬ 
meters);  discrete-layer  theory  (with  piecewise  linear  variation  of  the 
in  plane  displacements  in  the  thickness  direction);  simplified  discrete  layer 
theory  with  the  continuity  of  transverse  stresses  imposed  at  layer  interfaces 
to  reduce  the  number  of  displacement  parameters  to  five;  and  a  predictor- 
corrector  approach,  used  in  conjunction  with  the  first  order  shear- 
deformation  theory  (with  five  displacement  parameters  in  the  predictor 
phase). 

Based  on  the  numerical  studies  conducted,  llte  predictor-corrector 
approach  appears  to  be  the  most  effective  among  the  eight  modeling 
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approaches  considered.  For  antisymmetrically  laminated  rectangular 
plates  the  response  quantities  obtained  by  the  predictor  eorrector 
approach  ire  shown  to  be  in  close  agreement  with  exact  three-dimensional 
elasticity  solutions  for  a  wide  range  of  lamination  and  geometric  para¬ 
meters.  The  potential  of  this  approach  for  piedicting  the  response  of 
multilayered  anisotropic  plates  with  complicated  geometry  is  also 
discussed. 
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Stress,  Vibration,  and  Buckling 
of  Multilayered  Cylinders 

By  Ahmed  K.  Noor,1  Member,  ASCE,  and  Jeanne  M.  Peters1 

Abstract:  An  efficient  computational  procedure  is  presented  for  reducing  the 
cost  of  the  stress,  free  vibration,  and  buckling  analyses  of  multilayered  composite 
cylinders.  The  analytical  formulation  is  based  on  the  linear  three-dimensional  the¬ 
ory  of  elasticity  The  cylinders  are  assumed  to  have  simply  supported  curved  edges, 
and  the  fibers  of  the  different  layers  are  either  in  the  circumferential  or  longitudinal 
direction  The  fundamental  unknowns  consist  of  the  sis  stress  components  and  the 
three  displacement  components  of  the  cylinder.  Each  of  the  variables  is  expressed 
in  terms  of  a  double  Fourier  scries  in  the  longitudinal  and  circumferential  coor¬ 
dinates,  and  a  two  field  mixed  finite  element  model  is  used  for  the  discretization 
in  the  thickness  direction.  T  he  cylinder  response  associated  with  a  range  of  Fourier 
harmonics  in  the  longitudinal  and  circumferential  directions  is  approximated  by  a 
linear  combination  of  a  few  global  approximation  vectors,  which  are  generated  at 
panicular  values  of  the  Fourier  harmonics,  within  that  range.  The  full  equations 
of  the  finite  element  model  are  solved  for  only  a  single  pair  of  Fourier  harmonics, 
and  the  response  corresponding  to  the  other  Fourier  harmonics  Is  generated  using 
a  reduced  system  of  equations  with  considerably  fewer  degrees  of  freedom 
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Three-Dimensional  Solutions  for 
Antisymmetrically  Laminated 
Anisotropic  Plates 

Analytic  three-dimensional  elasticity  solutions  are  presented  for  the  stress  and  free 
vibration  problems  of  multilayered  anisotropic  plates.  The  plates  are  assumed  to 
have  rectangular  geometry  and  antisymmetric  lamination  with  respect  to  the  middle 
plane.  A  mixed  formulation  is  used  with  the  fundamental  unknowns  consisting  of 
the  six  stress  components  and  the  three  displacement  components  of  the  plate.  Each 
of  the  plate  variables  is  decomposed  into  symmetric  and  antisymmetric  components 
in  the  thickness  direction,  and  is  expressed  in  terms  of  a  double  Fourier  series  in  the 
Cartesian  surface  coordinates.  Extensive  numerical  results  are  presented  showing 
the  effects  of  variation  in  the  lamination  and  geometric  parameters  of  composite 
plates  on  the  importance  of  the  transverse  stress  and  strain  components. 
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A  K.  Noor  and  J.  M.  Print* 

George  Washinglon  University.  NASA  Langley  Research  Center,  Hampton.  VA  2J66J.  U  S. A. 

( RrcrlrrH  <y) 

Abstract— Two  recent  applications  of  (he  combined  use  of  mixed  formulation  with  splitting  methodolo¬ 
gies.  and  reduction  methods  are  presented.  Discussion  focuses  on  how  the  combination  can  significantly 
improve  the  efficiency  of  the  computational  process  and  enhance  the  physical  understanding  of  the 
response  Numerical  examples  are  presented  which  demonstrate  the  unique  features  and  potential  of  mixed 
formulations  in  advanced  analysis  systems. 
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BUCKLING  AND  POSTBUCKLING  ANALYSES  OF 
LAMINATED  ANISOTROPIC  STRUCTURES 


Ahmed  K.  Noor*  and  Jeanne  M.  Petersf 
George  Washington  University  Center,  MS-269 
NASA  Langley  Research  Center,  Hampton,  VA  23665,  U.S.A. 

ABSTRACT 

A  review  is  given  of  recent  advances  in  two  aspects  of  the  numerical  simulation  of  the  buckling 
and  postbuckling  responses  of  composite  structures.  The  first  aspect  is  exploiting  nontraditional 
symmetries  exhibited  by  composite  structures;  and  strategies  for  reducing  the  size  of  the  model 
and  the  cost  of  the  buckling  and  postbuckling  analyses  in  the  presence  of  symmetry-breaking 
conditions  (e  g.,  unsymmetry  of  the  material,  geometry,  and/or  loading).  The  second  aspect 
pertains  to  the  prediction  of  onset  of  local  delamination  in  the  postbuckling  range  and  accurate 
determination  of  transverse  shear  stresses  in  the  structure.  The  accuracy  and  effectiveness  of  the 
strategies  developed  are  demonstrated  by  means  of  a  numerical  example. 
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Novel  Computational  Strategies  for  Solution  of  Large-Scale 
Structural  Problems 

I 


Ahmed  K.  NOOR  and  Jeanne  M.  PKTIiRS 

George  Washington  University 
NASA  Langley  Research  Center 
I  lampion,  Virginia 


Summary 

Novel  computational  strategies  arc  presented  for  the  analysis  of  large  and  complex  structures. 
The  strategics  are  based  on  generating  the  response  of  the  complex  structure  using  large 
perturbations  from  that  of  a  simpler  model,  associated  with  a  simpler  structure  (or  a  simpler 
mathcmatical/diserete  model  of  the  original  structure).  Numerical  examples  are  presented  to 
demonstrate  the  effectiveness  of  the  strategies  developed. 
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ASSESSMENT  OP  COMPUTATIONAL  MODELS  FOR 
MULTILAYERED  COMPOSITE  CYLINDERS 

Ahmed  K.  Noor,  W.  Scull  Burton  and  Jeanne  M.  Pelers 
Oeorge  Washington  University 
NASA  Langley  Research  Center 
I  lampion,  Virginia  21665 


ABSTRACT 

A  study  is  made  or  the  effects  of  variation  in  Ihe  lamination  and  geometric  parameters  of 
multilayered  composite  cylinders  on  the  accuracy  of  the  static  and  vibrational  responses  pre¬ 
dicted  by  eight  modeling  approaches,  based  on  two-dimensional  shear  deformation  shell 
theories.  The  standard  of  comparison  is  taken  to  be  Ihe  exact  three  dimensional  elasticity  solu¬ 
tions,  and  ihe  quantities  compared  include  both  the  gross  response  characteristics  (e  g.,  vibration 
frequencies,  strain  energy  components,  average  through  the  thickness  displacements  and  rota¬ 
tions);  and  detailed,  duongli  (lie-thickness,  distributions  of  displacements,  stresses  and  strain 
energy  densities. 

Based  on  Ihe  numerical  studies  conducted,  a  predictor  corrector  approach,  used  in  conjunc¬ 
tion  with  Ihe  first-order  shear  deformation  theory  (with  five  displacement  parameters  in  Ihe 
predictor  phase),  appears  to  be  the  most  effective  among  Ihe  eight  modeling  approaches  con¬ 
sidered.  Eor  multilayered  orthotropic  cylinders  the  response  quantities  obtained  by  the 
predictor-corrector  appronch  are  shown  to  be  in  close  agreement  with  the  exact  three- 
dimensional  elasticity  solutions  for  a  wide  range  of  lamination  and  geometric  parameters.  The 
potential  of  this  approach  for  predicting  the  response  of  multilayered  shells  with  complicated 
geometry  is  also  discussed. 
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BUCKLING  AND  POSTBUCKLING  OF  T1  UN-WALLED 
COMPOSITE  STIFFENERS 

Ahmed  K.  Noor  and  Jeanne  M.  Peters 
*  George  Washington  University 

ABSTRACT 

A  review  is  given  of  the  different  types  of  symmetry  exhibited  by  the  prebnckling,  buck¬ 
ling  and  postbuckling  responses  of  composite  stiffeners.  A  simple  and  efficient  computational 
strategy  is  presented  for  generating  both  the  postbuckling  response  and  the  sensitivity  deriva¬ 
tives,  with  respect  to  lamination  and  material  parameters,  in  the  presence  of  symmetry-breaking 
conditions.  The  potential  of  the  foregoing  strategy  for  solving  practical  space  structures  prob¬ 
lems  is  discussed,  and  its  effectiveness  is  demonstrated  by  means  of  a  numerical  example  of  the 
postbuckling  and  sensitivity  analyses  of  a  composite  stiffener  with  a  zee-section  subjected  to 
uniform  end  shortening. 
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PRRDKTOR CORRP.CTOR  PROCP.DURP.S  FOR  STRRSS  AND 
RRRR  .VIBRATION  ANALYSPS  OP  MULTILAYPRKD 
^  COMPOSITE  PL  A  IPS  AND  SIIPLLS 

Ahmed  K.  Noor.  W.  Scott  Burton  nnd  leanne  M.  Peters 
(leorge  Washington  University 
NASA  Langley  Research  Cents*,, 

Hampton,  Virginia  2.1665 

ABSTRACT 

A  study  is  made  of  two  predictor-corrector  procedures  for  the  accurate  determination  or  the 

global,  as  well  as  detailed,  static  and  vibrational  response  characteristics  of  plates  and  shells. 

V 

Both  procedures  use  first-order  shear  deformation  theory  in  the  predictor  phase,  but  differ  in  the 
elements  of  the  computational  mode!  being  adjusted  in  the  corrector  phase.  The  first  procedure 
calculates  a  posteriori  estimates  of  the  composite  correction  factors  and  uses  them  to  adjust  the 
transverse  shear  stiffnesses  of  the  plate  (or  shell).  The  second  procedure  calculates  a  posteriori 
the  functional  dependence  of  the  displacement  components  cm  the  thickness  coordinate.  The 
corrected  quantities  are  then  used  in  conjunction  with  the  three-dimensional  equations  to  obtain 
belter  estimates  for  the  different  response  quantities.  Rxlensive  numerical  results  are  presented 
showing  the  effects  of  variation  in  the  geometric  and  lamination  parameters  for  antisymmetri¬ 
cally  laminated  anisotropic  plates,  and  simply  supported  multilayered  orthotropic  cylinders,  on 
the  accuracy  of  the  linear  static  and  free  vibrational  responses  obtained  by  the  predictor-corrector 
procedures.  Comparison  is  also  made  with  the  solutions  obtained  by  other  computational  models 
based  on  two  dimensional  shear  deformation  theories.  Por  each  problem  the  standard  <£  com¬ 
parison  is  taken  to  be  the  analytic  three-dimensional  elasticity  solution.  The  numerical  examples 
clearly  demonstrate  the  accuracy  and  effectiveness  of  the  predictor  corrector  procedures. 
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